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Abstract. We present observations of GRB 060124, the first event for which both the prompt and the afterglow emission could 
be observed simultaneously and in their entirety by the three Swift instruments. Indeed, Swift-BAT triggered on a precursor 
~ 570 s before the main burst peak, and this allowed Swift to repoint the narrow field instruments to the burst position ~ 350 s 
before the main burst occurred. GRB 060124 also triggered Konus-Wind, which observed the prompt emission in a harder 
gamma-ray band (up to 2 MeV). Thanks to these exceptional circumstances, the temporal and spectral properties of the prompt 
emission can be studied in the optical, X-ray and gamma-ray ranges. While the X-ray emission (0.2-10 keV) clearly tracks the 
gamma-ray burst, the optical component follows a different pattern, likely indicating a different origin, possibly the onset of 
external shocks. The prompt GRB spectrum shows significant spectral evolution, with both the peak energy and the spectral 
index varying. As observed in several long GRBs, significant lags are measured between the hard- and low-energy components, 
showing that this behaviour extends over 3 decades in energy. The GRB peaks are also much broader at soft energies. This is 
related to the temporal evolution of the spectrum, and can be accounted for by assuming that the electron spectral index softened 
with time. The burst energy {E\s„ ~ 5 x lO'^ erg) and average peak energy (jEp ~ 300 keV) make GRB 060124 consistent with 
the Amati relation. The X-ray afterglow is characterized by a decay which presents a break at ?i, ~ 10' s. 
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1. Introduction 

The Swift Gamma-Ray Burst Explorer dOehrels etal]l2004 
was successfully launched on 2004 Nov. 20. Its payload in- 
cludes one wide- field instrument, the ga mma-ray Burst Alert 
Telescope (BAT; 'Barthelmv et al."2005al 15-350 keV energy 
band), and two narrow-field instruments (NFIs), the X-Ray 
Telescope (XRT; iBurrows et alJ 120051 0.2-lOkeV) and the 
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Ultraviolet/Optical Telescope (UVOT; iRoming et aP l2005a[ 
1700-6500 A). The BAT detects the bursts, calculates their po- 
sitions to ^ 3' accuracy and triggers an autonomous slew of 
the observatory to point the two narrow-field instruments, typ- 
ically within 100 s from the burst onset. The XRT can provide 
S 5" positions, while the UVOT further refines the afterglow 
localization to ~ 0'.'5. 

GRB 060124 is the first event for which the three Swift in- 
struments have a clear detection of both the prompt and the af- 
terglow emission. Indeed, Swift-BAT triggered on a precursor 
on 2006-01-24 at 15:54:52 UT, ~ 570 s before the main burst 
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peak. This allowed Swift to immediately repoint the NFls and 
acquire a pointing towards the burst ~ 350 s before the burst oc- 
curred. The burst, which had a highly structured burst profile, 
comprises three major peaks following the precursor and had 
one of longest total durations (even excluding the precursor) 
recorded by either BATSE or Swift. 

GRB 060124 also triggered Konus-Wind dAptekar et alJ 

Il995h 559.4 s after the BAT trigger ( Golenets kiTet al.ll2006l) . 
The Konus light curve confirmed the presence of both the 
precursor and the three peaks of prompt emission. The main 
peak of GRB 060124 was also bright enough to trigger the 
FREGATE instrument aboard HETE-Il (HETE trigger 4012, 
ILamb et al. 2006) 557.7 s after the BAT trigger. 

The prompt emission of GRB 060124 was observed si- 
multaneously by XRT with exceptional signal-to-noise (S/N) 
and was detected by UVOT at V = 16.96 + 0.08 {T + 183 s) 
and V = 16.79 + 0.04 {T + 633 s). This fact makes it an 
exceptional test case to study prompt emission models, since 
this is the very first case that the burst proper could be ob- 
served with an X-ray CCD with high spatial resolution imaging 
down to 0.2 keV. There have been a handful of previous opti- 
cal detections while the burst prompt emission was still active, 
due to the increasing number of robotic foUowup telescopes . 
GRB 990123 was observed by ROTSE ( Aker lofet al] ll99^). 
and its light curve showed a flare which rapidly decayed, un- 
corrected with the gamma-ray emission. This behaviour was 
genera lly interpre ted as emission from the external reverse 
shock dSari & Piraa 1999b: Meszaros & Rees 1999). A reverse 
shock was also likely observed in GRB 050904 dBoer et alJ 
|2005). A different behaviour was observed for GRB 041219A 
iVesirand et al.„.2005; Blake et al., .20051) . In this case, the op- 
tical emission showed a temporal variabiUty correlated with 
the gamma-ray emission, suggesting a common origin for 
the two components . The prompt afterglow s of GRB 050319 
I Ouimbv et al Il2006t IWozniak et alJEool and GRB 050401 
I Rvkoff et al.l '200.5) showed a decay matching the late-time 
one, without the nee d for extra components. Last, GRB 050801 
jRvkoff et alJl2006l) showed a period of plateaux up to a few 
minutes after the burst, excluding contribution from the internal 
shocks, but possibly arguing for energisation from the central 
engine correlated with the X-ray emission. 

GRB 060124 was extensively followed from ground-based 
observatories with detections starting from one hour after 
the trigger The opti cal counterp art was identified by the 
Tautenburg telescope (Kann 2006). Further observations were 
can-i ed out with the ARIES 1.04m (Misra 2006). the HCT 
2m |Bhattetal.l [2006). the Cr AO 2.6m jRumvantsev et alJ 
I2OO6I) and the Asiago 1.82m (Masetti et al. 2006). Limits 
were given by ART (Torii 2006). MA STER (Lipunov et al 
| 2006l) . UoM 30cm (Sono da et alJl2006l). TARO T dKlotz et all 
12006 ^. and Loiano- 1.52m jGreco et all 12006). An absorp- 
tion redshift of z = 2.297 + 0.0 01 was determined based 
on the CIV/1^1548.1550 doublet^ ilMh^hal & Halperni l2006t 
ICenko et al."2006':'Prochaska et al."2006l). 

This paper is organized as follows. In Sect. |2] we describe 
our observations and data reduction; in Sect.|3]we describe our 
instrument-based spatial, timing and spectral data analysis; in 
Sect. 0] and Is] we describe our detailed, multi- wavelength tim- 



ing and spectral analysis of the prompt and afterglow emission. 
In Sect.|6lwe discuss our findings. Finally, in Sect. we sum- 
marise our findings and conclusions. Throughout this paper the 
quoted uncertainties are given at 90% confidence level for one 
interesting parameter (i.e., Ax^ - 2.71) unless otherwise stated. 
Times are referred to the BAT trigger Tq, t - T - To, unless 
otherwise specified. The decay and spectral indices are param- 
eterized as F(v,t) oc r"v^^, where Fy (erg cm"^ s"' Hz"') is 
the monochromatic flux as a function of time f and frequency 
v; we also use F = y6 H- 1 as the photon index, A^(£') « E^^ (ph 
keV ' cm"^ s '). We adopt a standard cosmology model with 
Ho = 70 km s-i Mpc"', Qm = 0.3, Qa = 0.7. 

2. Observations and data reduction 

2.1. BAT observations 

At 15:54:52 UT, Swift-BAT triggered and localiz ed GRB 
060124 (trigger 178 750, Holl and et al. 2006 : Fenimore et alJ 
2006). The on-board location was distributed 21 s after the 
trigger, and its position was RA(J2000) = 05''08"10% 
Dec(J2000) = -h69°42'33", with an uncertainty of 3'. The pai-- 
tial coding was 62%. The spacecraft started to slew to the burst 
position at f = 19 s, and settled at f = 89 s. Table n]reports the 
log of the Swift-BAT observations that were used for this work. 
Since BAT was triggered by the precursor, the event-by-event 
data which are available from To - 40 s < f < To -H 300 s only 
include the precursor and do not include the main burst. 

The temporal analysis around the main peak was per- 
formed using the mask-tagged 4 channel light curves (15- 
25, 25-50, 50-100, and 100-350keV bands) generated by the 
flight software. For the spectral analysis, the detector plane 
histogram (DPH) data (80-channel spectral data with time- 
intervals shown in the bottom panel of Fig.[3 were used. The 
BAT mask-weighted light curve (Fig. [0 shows a precursor 
from t - -3 to f = 13 s, with a peak count rate of ~ 1700 
counts s"' (15-350 keV) at the trigger time, then three major 
peaks: t - 520-550 s, t - 560-580 s (with the largest flux), and 
t = 690-710S. The inset of Fig. [T]shows the 16-s precursor. We 
note that this is the longest interval («; 500 s) record ed between 
a precursor and the main GRB event dT,a77atil2005l) . 

2.2. Konus observations 

The most intense part of GRB 060124 triggered Konus-Wind 
on 2006-01-24 at Tkw =16:04:13.9, 559.4s after the BAT trig- 
ger (taking into account the ~ 3-s propagation delay from Swift 
to Wind). It was detected by the S2 detector which observes 
the North ecliptic hemisphere; the incident angle was 43^4. 
The count rates are recorded by Konus in three energy ranges: 
Gl (18-70 keV), G2 (70-300 keV), and G3 (300-1 160 keV). 
The early (trigger) part of the time history is recorded from 
7'kw-0.512s to rKw+229.632s with time resolution ranging 
from 2 to 256 ms. The data before rKw-0.512s are collected 
in the waiting mode with 2.944 s time resolution. Time history 
recorded in the three energy ranges can be considered as 3- 
channel spectra. The full Konus time history of GRB 060124 is 
reported in Fig.|2l which shows a weak (4-cr and 2.5-<t in the 
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Fig. 1. Top: Mask-tagged (i.e., background subtracted) BAT 
light curve of GRB 060124 in the 15-350 keV energy band. 
The dot-long dash vertical line denotes the Konus trigger time 
Jkw- The dashed vertical lines denote the boundaries of the 
interval (A) over which we extracted the spectra in Sect. 13.11 
Inset: Mask-tagged BAT light curves during the precursor, ex- 
tracted in the standard 4 energy bands, 15-25, 25-50, 50-100, 
100-350keV, at a 3.2 s time resolution. The dot-short dash ver- 
tical line denotes the BAT trigger time Tq. Bottom: The BAT 
1 5-350 keV light curve around the time of the brightest part of 
the burst. The dashed vertical lines denote the boundaries of the 
five intervals over which we extracted the spectra in Sect. 13. II 

18-70keV and 70-300keV energy bands, respectively) pre- 
cursor corresponding to the Swift-BAT trigger time, no sta- 
tistically significant emission in any Konus energy band from 
? = 20 to 340 s, resuming emission at f = 340 s, the main 
double-peaked pulse in the interval t = 550-590 s, and weaker 
pulses up to f = 800 s. 

2.3. XRT observations 

The XRT observations of GRB 060124 started on 2006-01-24 
15:56:36 UT, only 104 s after the BAT trigger, and lasted until 
2006-02-23 23:32:57 UT. Table^reports the log of the Swift- 
XRT observations that were used for this work, which include 
a total net exposure time of ~ 867 s in Windowed Timing (WT) 
mode and ~ 247 ks in Photon Counting (PC) mode, spread over 
a ~31 d baseline. The monitoring is organized in 24 sequences. 
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Fig. 2. Konus Hght curves of GRB 060124 in the 18-70keV 
(a), 70-300 keV (b), 300-1 160 keV (c) and 18-1 160 keV (d) 
energy ranges at 2.944 s time resolution. The background lev- 
els are marked with dotted horizontal lines. The vertical lines 
denote the BAT (To, dot-short dash) and Konus (Tkw, dot-long 
dash) trigger times. 



The first (000) was performed as an automated target (AT) and 
with XRT in auto state. Therefore, in sequence 000 the XRT 
took an initial 2.5 s image (IM at f = 104 s) followed by WT 
frames (at t = 112 s) for the remainder of the first snapshot 
(continuous pointing at the target). During most of the first se- 
quence, the count rate was sufficiently high for XRT to remain 
in WT mode, the exception being 15 s recorded in PC mode 
(for a description of read-out modes, see iHill et alJl2004 . The 
latter were u sed to provide a first localization of the afterglow 
dMangano et al. 2006). The following observations were per- 
formed in both PC and WT modes, but the WT data have a 
very low S/N and were therefore not used here. 

2.4. UVOT observations 

Table ^ reports the log of the Swift-UVOT observations that 
were used for this work. The first UVOT observation was taken 
in IMAGE mode, because the burst interrupted a Target-of- 
Opportunity (ToO) pointing which had required this observing 
mode. Therefore, the UVOT did not execute the standard AT 
observing sequence, which uses EVENT mode to provide high 
time resolution. By contrast, the IMAGE mode data were col- 
lected in long integrations. After the first observation, the nor- 
mal follow-up mode was used. Figure|3shows the observations 
for which a detection was obtained (the upper limits are not 
shown because they do not usefully constrain the behaviour). 

3. Data analysis 

3.1. BAT and Konus analysis 

The BAT data were analysed using the standard BAT analysis 
software distributed within FTOOLS v6.0.4. The BAT ground- 
calculated position, based on data from t - -40 to -i-19s, is 
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Fig. 3. UVOT light curves. The V and B magnitudes are shown 
as (red) filled circles and (blue) empty squares, respectively. 



RA(J2000) = 05'>08'"30^6, Dec(J2000) = +69°43'27'.'9, with 
an uncertainty of 1'6 (radius, statistical and systematic eiTors). 
Mask-tagged BAT light curves were created in the standard 4 
energy bands, 15-25, 25-50, 50-100, 100-350 keV (Fig. 
and in the total band 15-350 keV (Fig.[0 at 3.2-s time resolu- 
tion. The BAT spectra were extracted over the interval which 
includes the precursor (A in Fig.^ t - ~l-5 to 13.5 s), and in 
the other 5 time intervals. Response matrices were generated 
with the task batdrmgen using the latest spectral redistribu- 
tion matrices. For our spectral fitting (XSPEC vl 1.3.2) we con- 
sidered the 14-150keV energy range and applied an energy- 
dependent systematic eiTor vector'. 

The Konus data were processed using standard Konus anal- 
ysis tools. From Tkw to Tkw -1-491.8 s, 64 spectra in 101 chan- 
nels were accumulated on timescales varying from 64 ms near 
the trigger (first four spectra) to 8.2 s afterwards. We also de- 
rived the Konus 3-channel spectra and the corresponding re- 
sponse for several time intervals prior to Tkw, as well as for 
the total burst, for the joint fitting with BAT and XRT. 

Tyo and T^o (the timescales over which 90% and 50% of the 
burst fluence is measured) were calculated from the BAT light 
curves with the task battblocks, using data binned to 1.6 s 
resolution in the time range t = 19.2-1 194.4s (thus excluding 
the precursor) and are reported in Table|2l These were also cal- 
culated for Konus data with Konus-specific tasks. We note, that 
the Konus T<)o and Tso in the G3 energy band are affected by 
low S/N effects and therefore are probably underestimated. 

Although the individual peak durations are comparable to 
the duration of other long GRBs, the total duration (excluding 
the pr ecursor) is among the longest recorded by either BATSE 
jPaciesas et al.1 19991) or Swift. In the 15-150 keV band the flu- 
ences of the precursor emission, the main emission and the fol- 
lowing peaks are (4.7 + 0.5) x lO"'', (1.25 + 0.03) x 10"^ and 
(1.05 + 0.1 1) X 10"^ erg cm"^, respectively. 



' http://heasarc.gsfc.nasa.gov/docs/swift/analysis/bat_digest.html 
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Fig. 4. UVOT (V), XRT, BAT and Konus light curves. The XRT 
light curves are corrected for pile-up. The count rates have been 
normalized to the peak of each light curve and offset vertically 
for clarity. 

3.2. XRT analysis 

The XRT data were first processed by the Swift Data Center 
at NASA/GSFC into Level 1 products (event hsts). Then they 
were further processed with the XRTDAS (vl.7.1) software 
package, written by the ASI Science Data Center (ASDC) and 
distributed within FTOOLS to produce the final cleaned event 
lists. In particular, we ran the task xrtpipeline (vO.9.9) ap- 
plying calibration and standard filtering and screening criteria. 
An on-board event threshold of ~0.2keV was applied to the 
central pixel of each event, which has been proven to reduce 
most of the background due to either the bright Earth limb or 
the CCD dark current (which depends on the CCD tempera- 
ture). For our analysis we selected XRT grades 0-12 and 0-2 
for PC and WT data, respect ively (according to Swift nomen- 
cl ature ; iBurrows et alJ20()5al) . 

The X-ray counterpart was detected at the position 
RA(J2000) = 05'>08'"26!03, Dec(J2000) = h-69°44'26'.'7, with 
an estimated uncertainty of 3'.'5. This position was determined 
using the xrtcentroid task (vO.2.7) on the sequence 002, 
which is not affected by pile-up, and it takes into account the 
correction for the misalignment between the telescope and the 
satellite optical axis. It is 63" from the refined BAT position 
(Sect.EU)- 



p. Romano: Prompt and afterglow emission of GRB 060124 



5 



- 


(a) 0.2-1 


keV 1 XRT - 

1 \ ' \ ' \ ' \ — ' 


- 


^ \ 

(b) 1-10 


keV 








r 


-H 1 


1 1 1 1 1 1 1 1 

1 \C) L — iu/ U.tC— L 










200 


400 600 800 1000 
Time since BAT trigger (s) 



Fig. 5. X-ray WT mode background-subtracted light curves of 
GRB 060124. (a): Soft band (S, 0.2-1 keV). (b): Hard band 
(H, 1-10 keV). (c): Ratio of hard to soft count rates. 

During sequence 000 the count rate of the burst was high 
enough to cause pile-up in the WT mode data. Therefore, to 
account for this effect, the WT data were extracted in a rectan- 
gular 40x20-pixel region with a 4x20-pixel region excluded 
from its centre. The size of the exclusion region was deter- 
mined following the procedure illustrated in the Appendix. To 
account for the background, WT events were also extracted 
within a rectangular box (40x20 pixels) far from background 
sources. WT events were extracted in three energy bands, 0.2- 
lOkeV (total), 0.2-1 keV (soft, S) and 1-lOkeV (hard, H). We 
dynamically subtracted their respective backgrounds. Figure|5] 
shows the soft and hard (background-subtracted and corrected 
for Point-Spread Function, PSF, losses) light curves, as well as 
the ratio H/S, with the BAT trigger as the origin of time. We 
also created WT light curves in the 1^ and 4-lOkeV bands. 
These are shown in Fig. |4] Tgo and Tgo were calculated from 
the XRT/WT light curves and are reported in Table|2 

During most of sequence 001 the PC mode data were piled- 
up, as well. Therefore, we extracted the source events in an 
annulus with a 35-pixel outer radius (~ 83") and a 3-pixel inner 
radius. These values were derived by comparing the observed 
and nominal PSF ( Vauahan et al. 2006). Given the decreasing 
brightness of the source, for the PC data collected in sequence 
002 and the following ones we extracted the spectra in circular 
regions of decreasing radius (35 to 10 pixels) to maximise the 
S/N. PC background data were also extracted in a source-free 
circular region (radius 42 pixels). 

Ancillary response files were generated with the task 
xrtmkarf within FTOOLS, and account for different extrac- 
tion regions and PSF corrections. We used the latest spectral 
redistribution matrices in the Cahbration Database (CALDB 
2.3) maintained by HEASARC. 

3.3. UVOT analysis 

The UVOT located the afterglow of GRB 060124 at 
RA(J2000) = 05'>08'"25!859, Dec(J2000) = h-69°44'27'.'41, 



with an internal accuracy of O'.'Ol and an estimated absolute 
astrometric uncertainty of 0'.'5. This is I'.'l from the XRT co- 
ordinates (Sect. lT^ . To avoid contamination from aV - 15.4 
source located 7" from GRB 060124, we performed aperture 
photometry using a circular aperture with a radius of 2" cen- 
tred on the optical afterglow. A sky annulus of width 5" and 
inner radius 17'.'5 was used since this annulus encloses a large 
sample of sky pixels, and it avoids the bright sources near the 
optical afterglow. Aperture corrections were performed to con- 
vert the 2" photometry aperture into the standard aperture radii 
used to define the photometric zero points. Six isolated stars 
were used to compute aperture coiTections for each exposure. 
The RMS scatter in the aperture corrections is typically 0.02 
mag. 

The instrumental magnitudes were transformed to Vega 
magnitudes using the photometric zero points in the Swift- 
UVOT Calibration Database (CALDB 2.3). Colour terms have 
not been applied, but preliminary calibrations suggest that 
they are negligible. The adopted photometric zero points were 
ZPv = 17.88 ± 0.09, ZPb = 19.16 + 0.12, ZPuvmi = 17.29 + 

0. 23, and ZPyvwi = 17.77 ± 0.02. Table|3i"eports the photom- 
etry. The UVOT fluxes reported in Table|3]are monochromatic 
fluxes computed at the central wavelength of the appropriate 
filter. 

Figure|3]shows the Swift-UVOT V and B Ught curves for all 
observations where a detection of the afterglow was achieved. 
The first two V-band points, at 183 and 633 s, differ from each 
other by 0.2 mag, which corresponds to a 2.3-cr variation. The 
UVOT light curve shows a remarkable amount of flickering 
(also see Sect. 14. 2> . We note that since these observations were 
performed in IMAGE mode, finer time resolution is not avail- 
able for these data. Therefore, the actual, possibly variable be- 
haviour of the optical light curve during the prompt phase re- 
mains unknown. We can estimate that if the prompt emission in 
the optical had occurred only in a time interval comparable to 
the one of the X-ray (and with no detectable emission outside, 

1. e. the effective exposure time is shorter), then the UVOT mea- 
surement would have been 1 .7 mag brighter. Even with the (ad- 
hoc) estimate for this shorter interval the optical flux would not 
be larger by the factor ~ 100 by which the X- and gamma-ray 
flux increase in the same amount of time. In order for the op- 
tical emission to fully track the hard energy band behaviour, a 
timescale of emission much shorter (£ 10 s) than that observed 
in the X-ray would have been required. 

4. Multi-wavelength timing 

4. 1 . The prompt phase 

The prompt phase of GRB 060124 was observed by UVOT, 
XRT, BAT and Konus, as shown in Fig.|?] which illustrates the 
prompt emission in different energy bands. A few observational 
features can be noticed. 

First, the relative importance of the three main peaks (as 
can be represented by the count rate ratios) varies through the 
prompt phase, in the sense that the third peak flux increases, 
relative to the second, as the energy band becomes softer. For 
example, the third peak in the 0.1-1 keV light curve is actually 
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Fig. 6. Shift of the emission peak with energy band. The XRT 
0.2-1 keV light curve is compared with the BAT 15-100 keV 
light curve. The lag between the peaks is ~ 10 s. 



stronger than the second one, which indicates that strong spec- 
tral evolution is taking plac e (also, se e Fig.l^];)- This has often 
been observed before (e.g., iFord et aljtl995^) . 

Second, the peak times of the emission shift with the energy 
band, as is shown in detail in Fig.|6l which shows the 0.2-1 keV 
XRT light curve compared with the BAT 15-lOOkeV one. In 
this extreme case, the lag between the peaks is ~ 10 s. 

The emission starts off simultaneously in all bands (albeit 
with the shift described above), although it has a considerably 
lower S/N in the BAT and Konus bands (see Fig.[2and Fig.|2] 
where the total light curves are shown). Furthermore, the emis- 
sion becomes less spiky as the energy of the band decreases. 
This is more quantitatively described by the XRT and BAT Tgo 
and Tso values which are reported in Tableland illustrated 
in Fig.0 While the Tgos are within 12% (fractional), the Tsqs 
vary by as much as 42%. This means that the peaks are much 
broader in the softer bands. 

The X-ray peaks can be formally modelled with 12 
Gaussians. When modelled with power laws, the prompt emis- 
sion peaks show very high rise slopes, ranging from 5 to 12 
if the time is referred to the BAT trigger, of the order of 1-2 
if the initial time is set to the first point of the flare. We note 
that there is evidence of continued flaring after the end of the 
first snapshot (see Figs.|5]and and our Tgos and T^os may 
therefore be underestimated. 



4.2. The afterglow phase 

The afterglow phase was observed by XRT and UVOT Since 
the XRT data are characterized by detection over a longer time 
and with a finer sampling, a joint fit is strongly biased in favour 
of the X-ray data. Therefore, we fit the XRT and UVOT data 
separately. 

The XRT/PC portion of the light curve (t > 10"^ s) was first 
fit with a simple power-law model using the BAT trigger as 
origin of time (we perform a fit with a different origin, T^w, 
below). The fit yields a decay index a = 1.36 + 0.02 and;^^^^^ = 




10 

Energy (keV) 



Fig. 7. Variation of Tgo (black filled circles) and T^o (red empty 
squares) with energy. The Tgo and T^q error bars correspond 
to the light curve binning (with the exception of the points at 
~ 700 keV). 



2.43 (94 degrees of freedom, d.o.f.). A fit with a broken power 
law F(t) = Kr"' for t < fb and F(t) = Kt^^"' (t/tbT"^ for 
f > fb, where fb is the time of the break, yields ai = 1.21 +0.04 
and a2 = 1.58 + 0.06, and a break at fb = (1.05+° |^) x 10^ s 
(X^^d - 1-87, 92 d.o.f.). This is a significant improvement over 
the simple power law (null hypothesis probability 1.9 x 10"^). 

We note that the XRT/PC light curve is affected by flick- 
ering, so that the residuals to the fits are dominated by rapid 
variations. Their intensity decreases as the time increases, a 
behaviour which has been observed before and which is due 
to both the logarithmic binning in our light curves and to the 
power-law decrease of the source counts. We tested the effect 
of flickering by removing the points in the light curve which 
were at least 2 cr above the model being tested for. Although 
the normalization of both the power-law and the broken power- 
law fits decrease, the slopes are consistent with the previous 
values within the uncertainties, and;^'^^^ drops down to ~ 1.2. 

It is also interesting to fit the light curves considering an 
initial time different from the BAT trigger time, i.e. the begin- 
ning of the main peak. To this end, we performed a fit of the 
PC data, using a broken power-law model with an initial time 
fixed to the Konus trigger time Tkw - 559.4 s. We obtain ai - 
1.20+°°^, a2 = 1.57+J}j;^, and break at fb = (1.05 +0.15) X lO^s 
(X^^^ = 1.88, 92 d.o.f.), which is consistent with what we found 
for the fits which adopted the BAT trigger as the starting time. 

The late-time UVOT light curve (f > 1 1 478 s) is character- 
ized by av = 0.72 ± 0.15 in the V band and as = 1.03 ± 0.23 
in the B band. A simultaneous fit of V and B band data with the 
constraint of a single slope yields obv = 0.82 + 0.06 (68% 
error, x^^^ - 1-13, 7 d.o.f.). We infer an observed value of 
the colour B - V = 0.47 + 0.08 and an optical spectral in- 
dex jSopt - 0.5 ± 0.3 (the latter was derived using an extinction 
correction of Ay = 0.447 mag). Figure |8] shows the homog- 
enized Swift-UVOT hght curve (the wavelength dependency 
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Fig. 8. XRT (filled circles) light curve in units of erg s~^ cm"^ 
and UVOT (empty squares) flux density (erg s"' cm"^ A'') 
scaled to match the X-ray. The solid line is the best fit model to 
the XRT afterglow data. The arrow marks the break position. 



of the UVOT bands was taken into account) superimposed (by 
scaling) on the XRT flux light curve. 

5. Broad-band spectroscopy: XRT, BAT and Konus 

5.1. The prompt phase 

As a way to characterize the fundamental burst parameters, 
we used the full Swift and Konus energy ranges. To this 
end we considered the 6 time intervals described in Sect. 13.11 
and extracted spectra from BAT in interval A, from XRT 
and BAT data in intervals 1 and 2, and XRT, BAT and 
Konus in intervals 3-5. We fit all the spectra extracted in 
the same time interval simultaneously. We considered an ab- 
sorbed power law model, with two absorption components, 

and one at the 



the Galactic one A^^ ,_— 

GRB redshift (z = 2.297. ICenko et al.ll2006l) . was free 
to vary between 7.65 x lO^" and 9.90 x 10^° cm^^ -phe 
lower limit was derived from the Galactic exti nction estimate 
along the GRB line of sight Av = 0.45 mag jSchlegel et alJ 
119981) . with a conversion in to Hydrogen column A^h = 1 .7 x 
lO^Uy cm-2 jPredehl & Schmitt 1995) . The second Umit is 
the up per limit of neutral Hydrogen maps jPickev & Lockmar] 
Il990 l). calculated with the nh tool within FTOOLS. We note 



that the Leideri/ Arge ntine/Bonn (LAB) Survey of Galactic HI 
jKalberla et alJ2005l) provides = 9.25 x 10^° cm ^ which 
is within the considered range. We considered different emis- 
sion models, ranging from a simple power-law, N(E) cc E^^, 
a cutof f power-law, A^(£ ') cc £'"'^ exp(-£'/£'o) and the GRBM 
model llBand et alJIl^gl . N(E) oc £-^1 exp(-£'/£'o) for E < 
(Fa - rO^o, N(E) oc £-^2 for £ > (Ta - fO^o, where £0 is 
the cutoff energy, and Fi and F2 are the low- and high-energy 
photon indices. The fits were performed in the 0.5-10, 14-150, 
20-2 000 keV energy ranges for XRT/WT, BAT and Konus, re- 
spectively. 

Our results for the best fitting model in each time interval 
are reported in Tableland illustrated in Fig.|9l At this time res- 



O.B 
0.8 
0.4 

0.3 


4 

3 

2 

1 
2 



1 

300 
200 
100 - 





(d) 



200 400 600 BOO 1000 

Time since BAT trigger (s) 

Fig. 9. Multi-wavelength time-resolved spectral analysis. The 
models and fit results are reported in Table 0] (a): BAT 15- 
350 keV light curve with the 6 time intervals chosen for the 
analysis marked; (b): in units of 10^^ cm"^; (c): photon 
index F; (d): peak energy £p. 



olution, the photon index generally increases (softens) through- 
out the prompt phase (Fig. |9j;). The precursor is significantly 
softer than the main burst prompt emission. However, we note 
that in the X-ray range, where a finer time sampling is available 
for spectroscopy, the hardness ratio tracks the count rate on a 
fine scale (see Fig.|5j;). Therefore, the emission is softer when 
its intensity is lower (also, see Sect. 15. 2> . There is marginal evi- 
dence for variations in the intrinsic column density, with a ~ 3- 
o" decrease between the third and fourth data point (Fig. |9j)). 
We note that the break energy lies outside the XRT band even 
in the last bin (Fig.|9|j). 

For comparison, we considered the Konus spectrum around 
the Konus trigger (f - 567.3-575.5) as representative of the 
most intense part of the burst. Fitting this spectrum with a GRBM 
model, we obtain Fi = l-13!o!7, ^2 = 2.12+JJ|^ ^nd £o = 
307![?o CrL = 0-84, 61 d.o.f.). We derive £p,Kw = 268!^^keV. 

We also fitted the total spectra (f = 301. 2-8 11. 2 s, last 
line in Table 0}. Based on this last fit, we derive an observed 
and rest-frame peak energy of £p = 193^3gkeV and £p - 
636^j2g keV, respectively. We also calculated integrated quan- 
tities, such as the fluence 'F, the peak photon and energy flux 
in several energy bands. These are reported in Table|5] 

The derived value of the isotropic energy, for a time in- 
tegration of 510 s in the source 1-10* keV energy range, is 
£iso - (4.2 + 0.5) x 10^^ erg. With these values, this bu rst is 
consistent with the £p-£iso relation f Amati et al1 l2002l) . The 
rest-frame isotropic luminosity in the 30-10000 keV range is 



-^p.iso 



(1.1 



with the £p-£p,iso 



0.1) X lO^^* erg s '. This value is also consistent 
relation CYonetoku et al.-2004) . 



5.2. Spectral evolution of the prompt phase with 
XRT/WT data 

To study the spectral evolution during the prompt phase in more 
detail than possible with BAT or Konus data, we created 23 
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Fig. 10. Time-resolved spectroscopy of the prompt emission 
of GRB 060124 with XRT/WT data, (a): Normalized XRT 
count rate, (b): Time evolution of the photon index. The ver- 
tical dashed lines mark the position of the peaks. 
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Fig. 11. Spectroscopy of the afterglow emission of 
GRB 060124 with XRT/PC data. Top: PC data fit with 
an absorbed power law model. The normalizations are relative 
to sequence 001, and are 0.29 (002), 0.09 (003), 0.05 (004), 
and 0.02 (005). Bottom: Ratio of the best-fit models to the 
data. 



time intervals that covered the prompt phase, and extracted 
XRT/WT spectra in the regions defined above. We fit them 
in the 0.5-lOkeV band, with an absorbed power law model, 
with the two absorption components described in Sect. l5.1l The 
results are shown in Fig. Qol As Fig. |5j; also shows, there is 
strong spectral evolution during the prompt phase. In particu- 
lar, the photon index has a minimum (i.e. the spectrum becomes 
harder) at the peaks of emission. Also, we note that the pho- 
ton index achieves a plateau following each peak in the light 
curve, with each plateau becoming successively softer. During 
this time interval, the photon index increases by about 1, reach- 
ing a value consistent with the later afterglow spectrum by the 
end of the last peak. 

5.3. XRT/PC spectrum of the afterglow 

The only spectroscopic information on the pure afterglow 
comes from the XRT/PC data. We performed a fit of the PC 
data in the 0.5-10 keV band, which were rebinned with a min- 
imum of 20 counts per energy bin to allow fitting within 
XSPEC. We considered an absorbed power law model, with 
the two absorption components described in Sect. 15.11 The 
fits for each sequence yield consistent photon indices and col- 
umn density at the GRB redshift, as reported in Table |5] We 
also fitted all data sets together and obtained a photon index 
rx,pc = 2.06 + 0.06, = 9.90 x 10^" cm^^, and a col- 
umn density at the GRB redshift A^^ - (LS+Jj ^) x 10^^ cm"^ 
(x'^^d ~ 1-07; 216 d.o.f., see Fig. II 1>. We confirm the presence 
of an excess of with respect to the Galactic Hydrogen col- 
umn along the line of sight on the order of 10^^ cm"^. 

5.4. Summary of the observed properties of 
GRB 060124 

1. The precursor/prompt phase. 



(a) For GRB 060124, we observed the longest interval (« 
500 s) recorded betw een a precursor and the main GRB 
event (lLazzatl2005l) . The precursor is spectrally softer 
than the following peaks (Table|4}. 

(b) The XRT, BAT and Konus light curves (0.2-2 000 keV) 
show the same overall structure consisting of three 
peaks. No time-resolved information is available from 
UVOT, but we observe that the optical behaviour during 
the prompt emission was significantly different from the 
X-ray one. 

(c) Based on the Tgos (x 300 s, Table|2j, although the indi- 
vidual peak durations are comparable to those of other 
long GRBs, the total durati on is among the lon gest 
recorded by either BATSE jPaciesas et alJ 1 19991) or 
Swift. 

(d) The peaks of the prompt emission shift with the energy 
band, with the peaks observed in the harder bands pre- 
ceding the ones observed in the softer bands. This lag 
is as much as ~ 10 s (between 15-100 keV and 0.2- 
1 keV). 

(e) Strong spectral evolution takes place during the prompt 
phase. We observe Ep moving from higher energies to 
lower energies (Sect. lsHI Fig. ll2t . The spectral evolu- 
tion follows a "tracking" behaviour (Fig.llOL with the 
spectrum being softer when the flux is higher The rela- 
tive importance of the three main peaks varies with the 
energy band; the third peak in the 0.1-1 keV light curve 
is actually stronger than the second one, as opposed to 
what is observed in all the other energy bands. 

(f) The peaks are much broader in the softer bands. 

(g) We derive a mean rest frame peak and isotropic energy 
E'p = 636^257 keV and Ei,o = (4.2 + 0.5) x lO" erg. 

2. The afterglow phase. 
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Fig. 12. Spectral energy distribution of GRB 060124 during the 
prompt phase. SEDs (1) and (2) are derived from simultaneous 
fits to XRT and BAT data fits, while SEDs (3)-(5) from XRT, 
BAT and Konus fits (see TablellJ. 

(a) The X-ray afterglow is modelled with a broken power- 
law decay with indices ai ~ 1 .2, 0-2^1-6 and a break 
at fb ~ 10^ s. 

(b) The optical afterglow is modelled with a power-law de- 
cay with index agy ~ 0.8 (the optical data are only 
available for t < 10^ s). 

6. Discussion 

One of the main results of Swift is the identification of a steep- 
decay phase which is observed in the early X-ray lig ht curves of 

Q GRBs, followed bv a shallower phase (.Tagliaferri et alJ 
iNousek etalJ 12004 lO'Rrien etal.l l200(il) . Ouite~ 
ten an erratic flaring behaviour is observed superimposed 
on these phases (^Burrows et al.l Eo05b: l alcone eLaL 2006: 
iRomano et al.l l2006l) . The most accredited explanation for 
the steep phase is that we are seeing high-latitude emis- 
sion due to the termination of the central engine activity 
jKumar & Panaitescul2000HZhang et alJ2006l) . while flares are 
probably due to the reactivation of the central engine. Here, for 
the first time, we have the chance to follow in the soft X-ray 
and in the optical/UV both the prompt phase and its transition 
to the afterglow phase. 

The panchromatic observation of GRB 060124 during the 
prompt phase provides the opportunity to study the broad- 
band emission of the GRB prompt emission in great detail. 
Time-resolved vFy spectra were analysed for 5 temporal seg- 
ments (Fig. I12> . each spanning nearly 4 decades in frequency 
(from the XRT to the BAT or Konus band). The flux lev- 
els at even lower frequencies (the UVOT band) are also mea- 
sured, but at a much coarser temporal resolution. Whenever the 
joint BAT-Konus data are available (bin 3-5), a peak energy 
Ep can be computed. The time-resolved spectral fits reveal an 
evolution of Ep to lower energies. The origin of has been 
extensively discussed in the literature (Meszaros et al. 1994; 
[Tavani. .1996: .Pilla & LoeK .1998: Zhang & Meszaro^. .2002; 



Pe'er&Waxmanl l2004 iRees & Mesz'irosI l2005t IPe'er etalJ 
20051: lRvdell2005l) . The suggestions range from the optically- 
thin characteristic synchrotron emission frequency, to the 
optically-thick Comptonised spectral peak, to the thermal peak 
of the fireball photosphere. In particular, there have been ques- 
tions about whether the synchrotron model can account for the 
observed spectrum ( iLlovd & Petrosianll2000l: iGhisellini et alJ 
l200Cl) . since the standard synchrotron model requires fast cool- 
ing so that the photon index below the Ep is expected to be 
~ 1 .5. This is because for internal emission models (e.g. the in- 
ternal shock model or the internal magnetic dissipation model), 
the magnetic fields in the dissipation region are so high that 
the electron cooling time scale is much shorter than the dy- 
namical time scale. The measured photon index in the bins 2-5 
are all close to this value, suggesting that the fast-cooling syn- 
chrotron model is a plausible mechanism at least for this burst. 
Lacking a temporally resolved UVOT light curve, one cannot 
derive a reliable time-resolved SED extending to the UVOT 
band. Nonetheless, connecting the time-averaged UVOT flux 
level to the flux level in the XRT band requires a further spec- 
tral break between the XRT and the UVOT band, which may be 
consistent with a cooling break or a synchrotron self-absorption 
break. 

Both the spectral peak energy Ep and the spectral indices 
evolve with time. This is likely due to the complicated in- 
terplay between particle heating and cooling as well as the 
emission contributions from many emission units (e.g. inter- 
nal shocks). One interesting property is that the second main 
peak (bin 5, Ep ~ 28 keV) is significantly softer than the 
first main peak (bin 4, Ep ~ 237 keV). Within the internal 
dissipation models, one h as Ep oc L'^^r"' oc L^^^y^^{6ty^ 
dZhang & MeszarosI l2002l) . where L is the wind luminosity 
of the GRB outflow, r is the emission radius, y is the bulk 
Lorentz factor, and 6t is the time variability from the central 
engine. The second part of proportionality is valid specifically 
for the internal shock model. The 15-150 keV fluence of the 
second main peak [(1.05 + 0.11) x 10"^ erg cm"^] is about 
one order of magnitude lower than that of the first main peak 
[(1.25 + 0.03) X 10"^ erg cm"^]. Given the similar durations of 
the two episodes, the luminosity of the second epoch is lower 
by a factor of ~ 10. This would account for a factor of 3 reduc- 
tion of the Zip. The data then require a larger emission radius 
(by a factor of ~ 3). In the internal shock scenario, this requires 
that the bulk Lorentz factor is somewhat higher. This may be 
consistent with the idea that the environment of the burst be- 
comes cleaner with time as the initial jet has loaded most of 
the baryons along the path. This would be the first piece of ev- 
idence of an increase of the Lorentz factor during the prompt 
phase. 

The fact that the pulses in the softer band tend to be broader 
is related to the temporal evolution of the spectrum. If the elec- 
tron spectral index does not evolve with time (as is usually as- 
sumed in the afterglow modeling), the rising/falling time scales 
of the light curves should be similar in different energy bands. 
However, if during the acceleration processes the electron spec- 
trum hardens with time, one can get a wider, soft wing with re- 
spect to a narrow, hard spike. The temporal decay slopes after 
the peaks are shallower than 2. This suggests that the decay is 
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not due to the adiabatic cooling of the emission region. Rather, 
electron heating may be still going on but with a reduced rate. If 
one assumes that the electron spectral index softens with time 
during the decay, one can account for a longer soft tail with 
respect to the narrow hard spike. This is supported by our find- 
ing that the photon indices are systematically smaller near the 
peaks (Fig. not . 

The super-long duration of the burst suggests that the 
GRB central engine is active for hundreds of seconds. This 
has become a general feature of both long and short GRBs 
in view of the commonly detected X-ray flares hundreds of 
seconds after the burst trigger, whic h are gener ally inter- 
preted as late central engine activities i Burrows et aL.2005 
Zhang et alJ f2006t iRomano et alJ 120061 iFalcone et al.l 1200^ 



Barthel mv et alJ l2005bt ICampana et all l2006t Liang et alJ 
2006) . The second main peak (excluding the precursor) would 



have been categorized as an X-ray flare had it been even 
softer The "precursor" emission leads the main peak by about 
500 s. This is inconsistent with precursor models that involve 
an erupting fireball from a massive stellar envelope, which 
predi cts a precursor leading the main episode by several sec- 
onds jRamirez-Ruiz et al .l2002h . A straightforward possibility 
is that the precursor is also due to early central engine activity. 
Suggestions to produce episodic, erratic cent ral engi ne activity 
include fragmentation of the collapsing s tar jKing et al..,2005 ) 
or the accretion disk jPerna et alJ l2006 t). a ma gnetic barrier 
modulated accretion flow (' Proga & Zhang ^'2006'). or the mag- 
netic activity of a massive neutron star central engine (Dai e t aL 
I2OO6I) . An interesting question for these models would be how 
to generate a weak precursor 570/(1 + z) = 173 s before the 
main episode. 

The origin of the optical emission during the prompt phase 
is uncertain from the available data. However, a tracking be- 
haviour would require timescales of optical emission in the 

order of ^ 10 s, wh ich is short, compared with GRB 990123 
J li^ — ' 

JAkerlof et al . 1999). Lacking a clear tracking behaviour with 
the prompt gamma-ray and X-ray light curves suggests that 
it may come from a dififerent emissio n site, e.g. the revers e 
shock ("Meszaros & Rees' 'l 997l Il999t ISari & PiranI |l999a'). 
Nonetheless, the level of the prompt optical emission clearly 
rules out the presence of a promine nt reverse shock com po- 
nent as discovered in GRB 990123 (Akerlofet al.'"l999') and 
GRB 021211 (Fox et al,.2003: Li et al. 20 03), which generally 
requires that the magnetic fields in the reverse shock region 
are sign ificantly s tronger than that in the for ward shock re- 
gion (iZhang et alJ >2003: Fan et al. 2002: Kumar & Panaitescul 
l2003l) . and that the fireball has comparabl e energy in baryonic 
mater ials and magnetic fields (i.e. cr ~ 1. 'Zhang & Kob avashil 
l2005h . The contribution of the reverse shock could be less 
significant when cr is low (baryonic) or ve ry high (Poynting- 
flux-dominated, Zhang & Kobayashi "2005V' In particular, var- 
ious physical processes could lead to the suppression of the 
reverse shock emission even if the fireball is baryonic (e.g 



lKobayas hi'2600: "Nakar & Piran''2004'; 'Kobayashi et al."2005) 
As a result, although the insignificant reverse shock emission 
is consistent with a Poynting flux dominated GRB model (e.g. 
iLvutikov & B landford' 2003: 'Roming et al.' '2005b"), the bary- 
onic GRB model is not disfavoured by the data. 



The early part of the afterglow (t < th) can be explained 
well in terms of the standard model dSari et alj [r998). In this 
case, the optical and X-ray data constrain the cooling frequency 
Vc between the optical and X-ray bands. In fact, the X-ray and 
optical spectral indices /3x - 1.1+0.10 and ySopt - 0.5 + 0.3 
diff'er by ^ 0.5, and constrain Vc ~ 2x 10'^ Hz. The temporal 
behaviour is fully consistent with this interpretation: both the 
optical and X-ray decay slopes are consistent with the model 
prediction: ax ~ 3/3x/2 - 1/2 = 1.15 + 0.10 (vaHd for v > Vc) 
and Qfopt ~ 3ySopt/2 = 0.75 + 0.4 (vaUd for v < Vc). The latter 
equation is valid only for a homogeneous external medium (the 
wind model would require a much steeper o-opt ~ 1.3). The 
electron distribution index is p - 2/3x = 2.2 + 0.15, which is in 
agreement with the shock theory. 

The X-ray decay slope after the break is too shallow for 
the break being due to a standard, sideways-expand ing jet, i n 
which case ax - p > 2 would be expected (Rhoads H 19991) . 
The break magnitude Aax = 0.39 + 0.08 may be consistent 
with a jet which propagates in a wind environment suff'ering no 
sideways expansion (for which Aa - 0.5). However, the wind 
model is ruled out by the early-time optical data. A solution 
which satisfies all the available constraints is a structured jet 
having a homogeneous core surr ounded by power -law wings 
with energy profile d/s/dQ oc ? (' Panaitescul2005l) . The break 
magnitude is dependent upon the index q (and on the location 
of Vc). With the observed values, we infer q ^ 0.85. The break 
magnitude below the cooling frequency (i.e. in the optical) is 
slightly different: Aofopt = 0.47. This difference, however, is 
quite small and may be easily missed. 

It is interesting to compare the expected break time ac- 
cording to the /iv.iso-^p-fb relation (Ghirlanda et al.. 20?)5t 
iLiang & Z hang 2005|, where th is the temporal break identified 
in optical data. According to this relation one would expect a 
break at fb = (2.1 + 1.2) x 10** s, for £^4^0 = (4.2 + 0.5) x 10^^ 
erg and E'^ = 636'^^^^ keV, which is broadly consistent with 
the observe d break time lO*" s. According to the Amati relation 
dAmati et al...2002.: .Amati .20061) . the expected intrinsic spec- 
tral break is £'p(Amati) = 583 keV. This is generally consistent 
with the measured value E'^ ~ 193!^^(1 + z) ~ 636^?^^ keV. 



7. Summary and conclusions 

We have presented observations of GRB 060124, the first event 
for which both the prompt and the afterglow emission could 
be observed simultaneously and in their entirety by the three 
Swift instruments. Since BAT triggered on a precursor ~ 570 s 
before the main burst peak, Swift repointed the narrow field in- 
struments to the burst position ~ 350 s before the main burst 
occurred. GRB 060124 also triggered Konus-Wind, which ob- 
served the prompt emission up to ~MeV energies. Thanks to 
these exceptional circumstances, we could report on the tem- 
poral and spectral properties of the prompt emission in the op- 
tical. X-ray and gamma-ray ranges, as well as on the optical 
and X-ray properties of the afterglow. 

We find that while the X-ray emission (0.2-10 keV) clearly 
tracks the gamma-ray burst, the optical component follows 
a different pattern, likely indicating a different origin, possi- 
bly the onset of external shocks. The prompt GRB spectrum 
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shows significant spectral evolution, with both the peak en- 
ergy and the spectral index varying. As observed in several 
long GRBs, significant lags are measured between the hard- 
and low-energy components, showing that this behaviour ex- 
tends over 3 decades in energy. The GRB peaks are shown to 
be much broader at soft energies. This is related to the tem- 
poral evolution of the spectrum, and can be accounted for by 
assuming the electron spectral index softened with time. The 
burst energy (Ei^o ~ 4.2 x 10^^ erg) and average rest-frame 
peak energy E'^ = 636+JgkeV) make GRB 010624 consistent 
with the Amati relation. The X-ray afterglow is characterized 
by a decay presenting a break at ~ 10^ s, possibly due to the 
crossing of the cooling frequency v^. 
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Table 2. Multi- wavelength and Tso- 



Obs/mode Energy range 


7*90 


T50 


(keV) 


(s) 


(s) 



XRT/WT 


0.2-1 


338+2 


160+2 


XRT/WT 


1^ 


302±1 


158±1 


XRT/WT 


4-10 


291±1 


138±1 


BAT 


15-25 


321±2 


60±2 


BAT 


25-50 


324±2 


54+2 


BAT 


50-100 


298±2 


50±2 


Konus 


300-1160 


47 ±9 


15 + 3 



Note: errors on and T^q are the time-resolution of the light curve used for the calculation. The Konus values are probably lower limits, 
due to low S/N effects. 



Table 3. UVOT multicolour data. 



Filter 


T(mid) 


T(expo) 


Magnitude^ 


Flux density'' 




(s) 


(s) 


(mag) 




V 


183.40 


153.31 


16.96 + 0.08 


0.522 ± 0.038 


V 


632.99 


725.22 


16.79 + 0.04 


0.610 + 0.022 


B 


1 1477.47 


139.43 


18.60 + 0.12 


0.220 + 0.024 


B 


17915.40 


884.89 


19.05 ± 0.06 


0.145 ±0.008 


B 


23674.52 


884.87 


19.51 ± 0.08 


0.095 ± 0.007 


B 


29314.24 


884.85 


19.53 ± 0.09 


0.093 ± 0.008 


V 


35202.82 


1573.78 


19.30 ± 0.09 


0.060 ± 0.005 


V 


41022.66 


1575.25 


19.44 + 0.10 


0.053 ± 0.005 


V 


46782.58 


1575.15 


19.52 + 0.10 


0.049 + 0.005 


V 


58332.60 


4725.17 


19.62 ± 0.07 


0.045 ± 0.003 


V 


75672.93 


4752.29 


19.86 ± 0.08 


0.036 ± 0.003 


V 


99028.58 


7947.60 


20.15 ± 0.08 


0.028 ± 0.002 



" Vega magnitudes. 

* In units of 10""' erg cm"^ s"' A"'. 

Note. Errors are 1-cr values. RMS errors in the photometric zero points have not been included. UVOT magnitudes were not corrected for 

Galactic extinction. Mid times are measured since the BAT trigger. 
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Time/Instrument 


Model A'jj Photon index 


£p 


Normalization 


xL (dof) 


Start 


End 


Selection" 


(10^2 cm-2) 


(keV) 


(ph cm-^ keV"') 


(s) 


(s) 



BAT Precursor 


PL 




1.8± 0.2 




(3.4±0.4) X 10-3(b) 


0.88 (57) 


-1.5 


13.5 


XRT-BAT Bin 1 


PL 


< 4.0 


1.17 ±0.07 




(7.2 ± 1.1) X 10-3(b) 


0.96 (70) 


301.2 


361.2 


XRT-BAT Bin 2 


PL 


T g+0.9 
^•^-0.8 


1.44 ±0.03 




(1.26 ± 0.09) X lO-^^(b) 


0.97 (137) 


361.2 


451.2 


XRT-BAT-Konus Bin 3 


CPL 




1.36 ±0.02 


J 32+29 


(3.0 ± 0.2) X 10-3(c) 


1.09 (404) 


451.2 


571.2 


XRT-BAT-Konus Bin 4 


CPL 


1.6 + 0.3 


1.45 ± 0.02 


9^7+153 


(1.7 ±0.1) X W-\c) 


0.99 (431) 


571.2 


691.2 


XRT-BAT-Konus Bin 5 


CPL 


1 9+0.3 


1 -70+0.03 


28+J° 


4.1!^ ^ X 10-* (c) 


0.83 (393) 


691.2 


811.2 


XRT-BAT-Konus (1-5) 


CPL 


1.9 ±0.2 


1.48± 0.02 




(1.19 ± 0.07) X IQ-Hc) 


1.06 (618) 


301.2 


811.2 



" SeeFig.|9| 

Normalization (K) at 50 keV. 
" Normalization (K) at 100 keV. 



Table 5. Multi-wavelength energetics of the prompt phase of GRB 060124. 



Energy band 


Energy fluence 


Peak photon flux 


Peak energy flux 


(keV) 


(erg cm"-) 


(ph cm-2 s-') 


(erg 


cm - s ) 


0.5-2 


(1.26 ±0.03) X 10-« 


15.8 ± 1.2 


(2.6 ± 


0.3) X 10-* 


2-30 


(7.3 ± 0.1) X 10-*^ 


19.1 ± 1.6 


(2.7 ± 


0.3) X 10-' 


30-400 


(1.8 ± 0.1) X 10-5 


7.8 ± 0.6 


(1.4 ± 


0.2) X 10-' 


25-50 


(3.60 ± 0.07) X 10-*^ 


3.2 ± 0.3 


(1.8 ± 


0.2) X 10-* 


50-100 


(4.7 ± 0.1) X 10^'^ 


2.5 ± 0.2 


(2.8 ± 


0.4) X 10-' 


100-300 


8.7+09 X 10-'' 


2.6 ± 0.2 


(7.1 ± 


0.9) X 10-' 


20-2000 


2.5:«-t X 10-5 


10.7 ± 0.9 


(2.2 ± 


0.3) X 10-" 


15-150 


(1.35 ± 0.03) X 10-5 


9.7 ± 0.8 


(7.8 ± 


1.0) X 10-' 


0.5-2000 


3.2+0^ X 10-5 


43.4 ± 3.5 


(2.4 ± 


0.3) X 10-« 



Table 6. Spectral fits of the afterglow phase of GRB 060124. 



Spectrum 


Model 


(1022 cm-2) 


Photon index 


Xl, (d.o.f.) 


Start time 
(s since Tq) 


End time 
(s tince Tq) 


Total exposure 
(s) 


XRT PC (001)" 
XRT PC (002)" 
XRT PC (003)" 
XRT PC (004)" 
XRT PC (mean AG)" 


PL 
PL 
PL 
PL 
PL 


1 6+" 

1 4+°-' 


2.14 ±0.13 
2.08 ± 0.09 
1.83 ±0.19 

2.06 ± 0.06 


1.14 (56) 
1.10 (104) 
0.82 (27) 
1.62 (10) 
1.07 (216) 


11408 

34398 

115809 

202206 

11408 


30726 

111665 

197945 

272587 

475924 


5479 

21979 

18495 

14269 

92360 



Model WABS*ZWABS (POWER) with a Galactic N° = 9.89 x 102° cm-2 ^j^^j j^gg photon index and column density at the GRB redshift. 
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Fig.A.l. Pile-up-driven spectral variations. XRT/WT photon 
index as a function of observed source count rate and size of the 
excluded region at the centre of the PSF (Sect.|Aj. The photon 
indices are shifted along the ordinate direction for clarity. 



Appendix A: Pile-up estimation in XRT/WT data 

The pile-up effect is a well known feature of CCD detectors 
llBalleil999l) . It occurs when the input count rate is so high that 
more than one X-ray photon is collected in a given 3x3 pixel 
region defining an X-ray event in a single frame and the charges 
produced by the two separate events are recorded as one. This 
effect causes groups of relatively soft X-ray photons to be de- 
tected as one single high-energy photon. The observed spec- 
trum of an astrophysical source, therefore, is distorted from its 
true form, with a general hardening of the slope (decrease of 
the photon index). 

Another feature of CCD detectors is that the charge pro- 
duced by a cosmic ray is usually spread over many pixels, as 
opposed to charge produced by X-ray photons. For this rea- 
son, a classification of the possible shapes of the charge dis- 
tribution is defined associating to each of them a grade num- 
ber: low grades correspond to confined charge distributions and 
high grades to extended charge distributions. Cosmic rays are 
then recognized and rejected on the basis of the event grade 
number Pile-up also occurs when two X-rays are collected in 
neighbouring pixels in one frame. In this case groups of low 
grade events can be recorded as a single high grade event. 
Then, when pile-up occurs, the observed grade distribution of 
the events will show a deficit of low grade events and an excess 
of high grade events as compared with event grade distribu- 
tion of fainter sources. Since higher grade events are rejected 
as cosmic rays, the pile-up effect causes an apparent loss of 
flux, particularly from the centre of the PSF where the count 
rate intensity is higher 

The standard way to correct for the pile-up effect is to de- 
termine the extension of the PSF core affected by pile-up and 
exclude from the analysis all the events that fell within that re- 
gion. Based on the three observable effects we have described, 
three different methods can be formulated to check for the pres- 
ence of pile-up and estimate the inner PSF region to exclude. 



They involve a study of the distortion of 1) the radial inten- 
sity profile; 2) the spectrum; 3) the grade distribut i on. Th e first 
method was discussed in detail by Vau^han et al.! ( 2006), with 
particular attention to Swift-XRT observations in PC mode, 
and it is widely applied as standard correction technique for 
PC data. The present analysis will focus on the following two 
methods as favourite pile-up diagnostic tools for Swift-XRT 
WT data. 

Ground testing and calibration of the XRT integrated sys- 
tem carried out in 2002 Sep. 23-Oct. 4 at the Panter laboratory 
of the Max-Planck-Institut fur extraterrestrische Physik showed 
that WT mode data are affected by pile-up for point source 
input intensities above 250 counts s"'. To determine the ex- 
tent to which pile-up distorts the spectral shape, we selected 5 
time intervals during which the observed source count rate was 
< 100, 100-200, 200-300, 300^00, and > 400 counts s'K For 
each count rate range, spectra were extracted from rectangular 
40x20-pixel regions (40 pixels long along the image strip and 
20 pixels wide) with a region of increasing size (0x20-15x20 
pixels) excluded from its centre. The spectra were fit with an 
absorbed power law with A^h fixed to the one obtained from 
the fit to the spectrum extracted with the exclusion of 15 pixels 
(which is assumed to be unaffected by piled-up). The spectra 
had 1433-2529 counts, depending on intensity cut and exclu- 
sion region. Figure IaTTI illustrates the pile-up-driven spectral 
changes. It demonstrates that in the range 0-100 counts s ' the 
source is not piled-up, since no point deviates from the horizon- 
tal line, which represents the best fit to the photon index values 
obtained in the 10-15 pixel exclusion range. For this intensity 
range, there is no need to exclude any pixels from the extrac- 
tion region. The source is only affected by moderate pile-up for 
100-300 counts s"'; in this case it is necessary to exclude 1 
pixel from the extraction region. Finally, two and four pixels 
must be excluded for 300^00 counts s"') and > 400 counts 
s"', respectively. 

To determine the extent to which pile-up distorts the grades, 
we studied the grade distribution for the same 5 time intervals 
and the same exclusion regions as for the spectral study. We de- 
rive analogous conclusions on the number of pixels to exclude 
from the extraction regions as a function of observed source 
count rate. The definitions of XRT grad es for timing niodes ar e 
different from the ones for PC mode ( Bur rows et alJ ^2005aV 
Grades in WT span from to 15, and good events are likely 
associated to grades from to 2 only. If we assume that a good 
approximation of the expected grade distribution with no pile- 
up is obtained excluding 10 central pixels, the distributions we 
obtained show a deficit at grade zero (associated to an increase 
at higher grades) that becomes negligible if a region of at least 
4-pixels is excluded. 
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Table 1. Observation log of GRB 060124. 



Sequence 


Obs/mode 


Start time (UT) 


End time (UT) 


Net exposure 


Time since trigger 






(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 


(s) 


(s) 


00178750000 


BAT/event-by-event 


2006-01-24 15:54:13 


2006-01-24 15:59:55 


342.3 


-40 


00178750000 


BAT/mask-weighted LC 


2006-01-24 15:55:12 


2006-01-24 16:14:51 


1179.2 


19.2 


00178750000 


BAT/DPH 


2006-01-24 15:59:54 


2006-01-24 16:08:24 


510 


301.2 


00178750000 


XRT/IM 


2006-01-24 15:56:36 


2006-01-24 15:56:38 


2.5 


104 


00178750000 


XRT/WT 


2006-01-24 15:56:44 


2006-01-24 16:11:34 


867 


112 


00178750000 


XRT/PC 


2006-01-24 15:58:13 


2006-01-24 15:59:01 


15 


200 


00178750001 


XRT/PC 


2006-01-24 19:05:00 


2006-01-25 00:26:58 


5501 


11408 


00178750002 


XRT/PC 


2006-01-25 01:28:10 


2006-01-25 22:55:56 


22067 


34399 


00178750003 


XRT/PC 


2006-01-26 00:05:00 


2006-01-26 22:53:56 


18569 


115809 


00178750004 


XRT/PC 


2006-01-27 00:04:58 


2006-01-27 19:37:58 


14327 


202206 


00178750005 


XRT/PC 


2006-01-28 00:08:40 


2006-01-30 23:17:57 


47631 


288829 


00178750007 


XRT/PC 


2006-01-31 00:24:14 


2006-02-01 23:24:57 


44595 


548963 


00178750009 


XRT/PC 


2006-02-02 00:41:55 


2006-02-02 08:42:03 


4759 


722824 


00178750010 


XRT/PC 


2006-02-03 23:24:51 


2006-02-03 23:35:58 


667 


891000 


00178750011 


XRT/PC 


2006-02-04 11:58:22 


2006-02-04 20:28:57 


4816 


936210 


00178750012 


XRT/PC 


2006-02-05 13:40:08 


2006-02-06 23:59:57 


10779 


1028716 


00178750013 


XRT/PC 


2006-02-07 01:11:20 


2006-02-07 23:49:58 


9483 


1156588 


00178750014 


XRT/PC 


2006-02-08 01:12:24 


2006-02-08 23:51:55 


10310 


1243053 


00178750015 


XRT/PC 


2006-02-09 06:07:56 


2006-02-09 23:56:57 


6443 


1347185 


00178750016 


XRT/PC 


2006-02-10 00:04:58 


2006-02-10 22:43:49 


423 


1411807 


00178750017 


XRT/PC 


2006-02-11 01:42:30 


2006-02-13 23:03:56 


7248 


1504059 


00178750018 


XRT/PC 


2006-02-15 00:10:07 


2006-02-15 22:55:56 


11664 


1844115 


00178750019 


XRT/PC 


2006-02-17 00:20:47 


2006-02-17 23:07:57 


15741 


2017556 


00178750000 


UVOT/V 


2006-01-24 15:56:38 


2006-01-24 15:59:14 


153 


106 


00178750000 


UVOT/V 


2006-01-24 15:59:17 


2006-01-24 16:11:34 


725 


265 


00178750001 


UVOT/B 


2006-01-24 19:04:58 


2006-01-24 19:07:20 


139 


11407 


00178750001 


UVOT/UVM2 


2006-01-24 19:35:31 


2006-01-24 19:39:15 


221 


13239^ 


00178750001 


UV0T/UVM2 


2006-01-24 19:39:18 


2006-01-24 19:40:03 


43 


13466^ 


00178750001 


UVOT/B 


2006-01-24 20:45:53 


2006-01-24 21:01:02 


885 


17461 


00178750001 


UV0T/UVW2 


2006-01-24 21:01:08 


2006-01-24 21:04:53 


221 


18376^ 


00178750001 


UVOT/UVW2 


2006-01-24 21:04:56 


2006-01-24 21:08:41 


221 


18604^ 


00178750001 


UVOT/UVW2 


2006-01-24 21:08:44 


2006-01-24 21:12:29 


221 


18832'' 


00178750001 


UVOT/UVW2 


2006-01-24 21:12:32 


2006-01-24 21:16:17 


221 


19060'' 


00178750001 


UVOT/B 


2006-01-24 22:21:52 


2006-01-24 22:37:01 


885 


23220 


00178750001 


UV0T/UVW2 


2006-01-24 22:37:07 


2006-01-24 22:40:52 


221 


24135'' 


00178750001 


UVOT/UVW2 


2006-01-24 22:40:55 


2006-01-24 22:44:40 


221 


24363^ 


00178750001 


UVOT/UVW2 


2006-01-24 22:44:43 


2006-01-24 22:48:28 


221 


24591^ 


00178750001 


UVOT/UVW2 


2006-01-24 22:48:31 


2006-01-24 22:52:16 


221 


24819^ 


00178750001 


UVOT/B 


2006-01-24 23:55:52 


2006-01-25 00:11:00 


885 


28860 


00178750001 


UV0T/UVW2 


2006-01-25 00:11:07 


2006-01-25 00:14:52 


221 


29775" 


00178750001 


UVOT/UVW2 


2006-01-25 00:14:55 


2006-01-25 00:18:40 


221 


30003" 


00178750001 


UVOT/UVW2 


2006-01-25 00:18:43 


2006-01-25 00:22:28 


221 


3023 P 


00178750001 


UV0T/UVW2 


2006-01-25 00:22:31 


2006-01-25 00:26:16 


221 


30459" 


00178750002 


UVOT/V 


2006-01-25 01:28:09 


2006-01-25 01:55:00 


1574 


34398 


00178750002 


UVOT/V 


2006-01-25 03:05:09 


2006-01-25 03:32:00 


1575 


40217 


00178750002 


UVOT/V 


2006-01-25 04:41:09 


2006-01-25 05:08:00 


1575 


45977 


00178750002 


UVOT/V 


2006-01-25 06:17:09 


2006-01-25 09:57:00 


4725 


51737 


00178750002 


UVOT/V 


2006-01-25 11:06:09 


2006-01-25 14:46:00 


4752 


69078 


00178750002 


UVOT/V 


2006-01-25 15:54:41 


2006-01-25 22:56:00 


7948 


86389 



Upper limit on detection. 



